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Semiconducting graphene from highly ordered substrate interactions 
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While numerous methods have been proposed to produce semiconducting graphene, 
a significant bandgap has never been demonstrated. The reason is that, regardless of 
the theoretical gap formation mechanism, disorder at the sub-nanometer scale prevents 
the required chiral symmetry breaking necessary to open a bandgap in graphene. In 
this work, we show for the first time that a 2D semiconducting graphene film can be 
made by epitaxial growth. Using improved growth methods, we show by direct band 
measurements that a bandgap greater than 0.5 eV can be produced in the first graphene 
layer grown on the SiC(OOOl) surface. This work demonstrates that order, a property 
that remains lacking in other graphene systems, is key to producing electronically 
viable semiconducting graphene. 


It is well known that the first graphene layer grown on 
the SiC(0001) surface is not electronic graphene. That 
is, the first “buffer” graphene layer does not show the 
linear dispersing 7r-bands (Dirac cone) expected at the 
iC-point of metallic graphene. m The lack of 7r-bands 
in experimental band maps of the buffer layer [2| sup¬ 
ported the theoretical conclusion that sufficiently strong 
covalent bonds between the buffer layer and the SiC in¬ 
terface would push the graphene 7r-bands below the SiC 
valance band maximum. w Aside from these very early 
studies, research on the SiC graphene buffer layer faded 
and was subsequently eclipsed by a wide variety of other 
unsuccessful ideas to open a band gap in exfoliated or 
chemical vapor deposition (CVD)-grown graphene. 0 

One method to open a band gap in graphene is by pe¬ 
riodic bonding to either all A or all B sites, which breaks 
graphene’s chiral symmetry (referred to as graphene 
functionalization). The buffer graphene, commensu- 
rately bonded to the SiC(0001) surface, should have been 
an excellent example of a functionalized system that in¬ 
duces a bandgap. Despite the buffer graphene’s potential 
to be functionalized by a commensurate and, most im¬ 
portantly, ordered array of Si or C atoms in the SiC, there 
was a major research shift to functionalize CVD-grown 
graphene. Efforts to functionalize CVD graphene by a 
number of other methods have been a major research 
area. As of this writing, no functionalized graphene, or 
graphene modified by any other proposed method, has 
been developed that produces a workable semiconducting 
form of graphene. The problem with these methods is the 
inherent disorder introduced in the functionalization [7] 
and growth process. [8] In fact, the lack of a graphene 
bandgap was the motivation to shift research to metal 
dichalcogenides despite the inability to grow them at the 
level of purity and order required for industrial scale elec¬ 
tronics. 

In this work, we use furnace-grown graphene to pro- 



k y ( A' 1 ) K(k-') 


FIG. 1. (a) An ARPES cut through the graphene K-point of 
an under-grown 6 a/ 3 surface. k y is perpendicular to T — K. 
The states g± and g<i (dashed lines) observed by Emtsev et al. 
[2] are marked, (b) The same cut as in (a) for growth 20°C 
higher. Circles mark the peak positions along part of the e\ 
band. 


duce a structurally well ordered buffer graphene on 
the SiC (0001) surface. Angle resolved photoemission 
(ARPES) measurements show new dispersing 7r-bands 
that are not observed in samples grown by previous meth¬ 
ods. These bands live above the SiC valance band max¬ 
imum near the Fermi Energy, Ep. The new band struc¬ 
ture is a result of improved order caused by a higher 
growth temperature which, for the first time, gives rise 
to a well ordered 6x6 reconstruction in surface x-ray scat¬ 
tering experiments. 9 The bandgap, which is >0.5 eV, 
appears to be the result of the chiral symmetry breaking 
caused by the 6x6 reconstruction. The ARPES show that 
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FIG. 2. A constant E—Ep = — 1 eV cut through part of a ML 
graphene BZ showing replica cones. Blue dots mark single 
6x6 umklapp scattering replicas of the Dirac cones (si,S 2 ). 
Umklapp scattering of the Dirac cones from SiC lxl (purple 
dots) and SiC lxl plus 6x6 reciprocal lattice vectors (green 
dots) are also shown. 


the buffer graphene layer on SiC is a true semiconductor, 
the goal of the first graphene electronics research. mm 

The substrates used in these studies were n-doped 
n = 2 x 10 18 cm -2 CMP polished on-axis 4H-SiC(0001). 
The graphene was grown in a controlled silicon sublima¬ 
tion furnace. 12 Graphene growth is a function of tem¬ 
perature, time, and crucible geometry that sets the sil¬ 
icon vapor pressure. With the current crucible design, 
a monolayer (ML) graphene film will grow at 1520°C in 
20 min. Using the same crucible, the semiconducting 
buffer layer discussed in this paper will grow at a tem¬ 
perature 160° C lower than the ML in the same amount of 
time. Growing 20°C lower than the optimum buffer tem¬ 
perature will still give the same (OV^xbv^RBO 0 LEED 
pattern (referred to as 6>/3 in the subsequent discussion) 
as the optimum buffer film but will not show the gapped 
7r-bands discussed below. 

Early ARPES work on the UHV-grown recon¬ 
structed SiC(0001) surface (referred to as the graphene 
buffer layer in later literature) found that two non¬ 
dispersing states gi and g 2 at -0.5 and -1.6 eV were the 
only band features between Ep and the SiC valance band 
maximum. 2 These states were interpreted as localized 
Mott-Hubbard states hybridized from SiC surface dan¬ 
gling bonds. We can reproduce these states by heating 
the SiC 20° C cooler than the optimal buffer growth tem¬ 
perature. Figure [lja) shows an ARPES cut through the 
graphene iC-point from this “sub-buffer” film. The two 
surface states seen in previous work are clearly visible. 
However, by heating 20° C higher, a new dispersing band, 
ei(fc) appears [see Fig. EF)]- The new surface state is ro¬ 
bust, being reproducible in multiple samples. Note that 
a faint linear Dirac cones appears at k y = 0. This is due 
to a small amount of ML graphene (<2%) that typically 


nucleates at intrinsic step edges. m The Dirac point of 
the partial monolayer has the typical n-doping (0.55eV 
below Ep). 

Another indication of the improved sample order is 
the quality of the monolayer grown above the optimum 
buffer. Figure [2] shows a constant energy cut through 
part of the Brillouin Zone (BZ) of a ML graphene film. 
In addition to the Dirac cone, replicas of the Dirac cone 
from umklapp scattering processes are also visible. All 
replica cones from the K th K -point can be indexed us¬ 
ing reciprocal lattice vectors of the SiC 6x6 unit cell: 
Gx(m,n) = msi ±ns 2 , where |si| = |s 2 I ~ 6 l a Sicl [see 
Fig. !• In the ordered ML films, replica cones are clearly 
seen from both l st -order in the 6x6 unit cell (si,s 2 ) 
and from multiple scattering processes involving 1 st -order 
(si,s 2 ) plus a SiC G vector (e.g., the Gx(7,0) and 
G^(7,l)). Early UHV grown samples only showed 1 st - 
order replicas (i.e, n,m=l).[14j The fact that so many 
ARPES replicas bands are observed in these films, along 
with the 6th order x-ray diffraction rods,[9j testifies to 
the film’s improved order. 

Detailed ARPES measurements from these improved 
samples show that the e\{k) band [Fig. [Tj(b)] is a gapped 
graphene 7r-band. Figure |3|a) shows a constant energy 
cut though part of the BZ of a buffer layer graphene near 
the ei(fc) band maximum. Three lobes are visible. These 
lobes represent a second dispersing bands, e 2 (fc), that is 
marked in the YKM' cut in Fig. [3](b) and (c). Again, a 
Dirac cone from a small amount of ML graphene is visi¬ 
ble. The two bands are independent of the perpendicular 
momentum k±(E) and therefore cannot be due to bulk 
bands. The tops of both bands lie A E ~ 0.5 eV below 
Ep , or 1.8 eV above the valance band maximum of SiC 
interface, indicating that the buffer is a wide band gap 
semiconducting form of graphene. A schematic of the 
two bands is shown in Fig. [4j The ei(fc) band appears as 
a gapped 7r-band that disperses slower perpendicular to 
YK than along either YK or KM directions [see Table [I]. 
The linear part of the ei(fc) band has a velocity, u, that is 
significantly lower than the Fermi velocity, vp, reducing 
to nearly half vp perpendicular to YK [see Table [T]. 


TABLE I. The band velocity (v) and effective mass (m*) near 
the K -point near the 7r-band maximum, m* is estimated 
assuming parabolic bands near Ef. 


Band 

v/vf 

m* /m e 

ML Dirac cone 

1.0 

- 

ci (± TK) 

0.55 ±0.01 

1.0 ±0.02 

ci {TK) 

0.63 ±0.1 

1.5 ±0.5 

ei (KM) 

0.80 ±0.1 

0.55 ±0.05 

e 2 (± TK) 

0.98 ±0.07 

0.25 ±0.02 

e 2 (TK) 

- 

1.5 ±0.1 


The e 2 (fc) band is 3-fold symmetric, extending towards 
T and dispersing perpendicular to YK. Figure[3|c) shows 
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FIG. 3. (a) A constant energy cut through the graphene BZ near the K-point (E — Ep — —0.41eV, kv — 70eV). Dashed lines 
mark the boundary of the BZ. (b) A cut through the surface bands in the TKM' direction. Circles mark the peak positions 
along part of the ei and e 2 band along with a few higher binding energy bands. A weak Dirac cone from a partial ML is shown, 
(c) (b) A cut perpendicular to TK through the 62 band [vertical black dashed line in (a)]. Circles mark the peak positions of 
the 62 band. 


a cut perpendicular through the lobe in Fig. [3^a) . The 
band velocity of 62 (h) perpendicular to TK is nearly the 
same as monolayer graphene [see Table [I]. The 6 \ band 
has an effective mass (m*) that ranges between 0.55 to 
1.5ra e , while 62 is a light band perpendicular to TK. 

In a broad sense, the gapped band structure near the 
iC-point strongly suggests chiral symmetry breaking that 
mixes the 7 r bands from the K and K' points. US] Any 
periodic potentials that break the AB symmetry in the 
graphene through bond formation, chemical or strain 
fields, or finite size effects can open a gap in graphene. 
Weak interactions like those in bilayer graphene only 
produce small gaps. m The strain necessary to induce 
Kekule distortions, which produce band gaps of the 
order observed in these samples, are large enough to 
tear the graphene. IS] That strain level is inconsistent 
with the 0.4% strain measured by X-ray scattering. mm 
Strong bonding using Aryl functionalization [7, 20] has 
given large gaps, but the disorder inherent during the 
functional group’s incorporation into the graphene lat¬ 
tice leads to poorly resolved band structure and low 
mobilities. [21] 

A theoretical understanding of the buffer layer, and 
therefore an understanding of the origin of the observed 
gap, is difficult because of the excessive calculation time 
associated with exploring different models for the large 
6^/3 unit cell. Rather than calculating the full 6>/3 cell, 
early calculations instead used an almost commensurate 
V3 x R30 SiC reconstruction to make calculations 
more tractable. mm These calculations predicted that 
strong sp 3 bonds between 2/3 of the interfacial Si atoms 


FIG. 4. A schematic representation of the a (k) and 62 (k) 
buffer layer bands near the top of the 7r-bands around the 
K- point. 


and the buffer buffer graphene caused the 7r-bands to 
shift above and below the conduction band minimum 
and valance band maximum, respectively. The calcula¬ 
tions also predicted a metallic, slightly delocalized, sur¬ 
face state near Ep due to the remaining unbounded SiC 
Si atoms in the interface, similar to the states observed 
experimentally in the earlier, less ordered samples like 
in Fig. Ba).E] These approximate models are clearly in¬ 
sufficient to explain the observed bands. Only one ab 
inito calculation by Kim et al.[22] has calculated the 
band structure for the buffer using a full 6\/3 cell. While 
the calculation was restricted to a bulk terminated SiC 
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FIG. 5. (a) The calculated graphene buffer layer (6\/3 X 

6^3)R30° cell on a relaxed bulk terminated SiC(OOOl) 
surface. [22] Green open circles are buffer carbon atoms that 
are bonded to the SiC surface. Black circles are carbon not 
bonded to substrate. Chains of atoms define superhexgon re¬ 
gions. (b) An outline of the calculated bands (red dashed 
lines) from a bulk terminated SiC plus buffer layer from 
Ref. [22] overlaid onto the experimental bands. The theo¬ 
retical bands have been shifted 0.13 eV lower to better match 
the 62 band. 


interface, [22 ^ it does give some insight into the origin of 
the observed gap when compared to the ARPES results. 

Kim et al. [22] find that about 25% of the carbon 
atoms in the buffer graphene are covalently bonded to 
Si atoms on the SiC interface. The resulting structure 
is a hexagonal network of graphene ribbon-like struc¬ 
tures with the remaining buffer carbon atoms covalently 
bonded to the SiC surface [see Fig. [5|a)]. Similar hexag¬ 
onal networks, either structural or electronic, have been 
observed in Scanning Tunneling Microscopy or produced 
theoretically. [23H25] Kim et al.’s DFT calculations show 
that the 7r-orbitals of carbon atoms on the superhexago- 
nal boundaries (or ribbons) give rise to several bands near 
the K -point above and below Ep. These bands are over¬ 
laid on our measured bands in Fig.[5][b). We have shifted 
the calculated bands by -0.13eV to match the 62 (k) band 
maximum. Like the experimental e\(k) and 62 (k) bands, 
the theoretical model shows that the covalent bonding 
to the SiC does not completely destroy the 7r-bands as 
earlier calculations predicted. Nonetheless, the calcula¬ 


tions do not reproduce several important features of the 
experimental bands. The calculations to not predict the 
formation of a band gap. They also do not correctly re¬ 
produce the dispersion of the e\ band, especially from T 
to K. 

We suggest that while the ribbon structure produced 
in the model of Kim et al.[22] may be correct, the large 
amount of covalent bonds associated with a bulk ter¬ 
minated SiC likely over estimates the graphene-SiC in¬ 
teraction. It is more likely that the buffer graphene is 
bonded to the SiC through a smaller number of sites. 
A lower number of graphene-Si bonds is more consis¬ 
tent with STM measurements that suggest the buffer 
lies above a small set of Si-t rimers. [26, A reduced buffer- 
SiC bonding geometry is also consistent with both x-ray 
scattering (271 and x-ray standing wave experiments, [28] 
which find a reduced Si-concentration and an increased 
C-concentration in the SiC layer below the buffer. We 
suggest that reduced substrate bonding would still be 
sufficient to strain the buffer to produce the observed rib¬ 
bon network. The ribbon network, which is now isolated 
from the substrate, would result in a finite size induced 
band gap, as observed. [29] 

In this work we show that, despite claims to the con¬ 
trary, a semiconducting form of graphene can be manu¬ 
factured. Using improved growth methods, we have been 
able to produce a well ordered single layer of graphene on 
the SiC(0001) surface. The first graphene layer, known as 
the buffer layer, is a semiconducting form of graphene. 
Using ARPES, we show that the improved sample or¬ 
der leads to new bands with a band maximum that lies 
AF^O.5 eV below the Fermi Level. Depending on where 
the conduction band lies, the bandgap of this form of 
graphene must be >0.5 eV. While no theoretical model 
predicts the measured bands, the experimental bands re¬ 
semble those from a network of graphene ribbons that 
are distortions in the buffer layer. The distortions would 
be due to strain strain relief in the film caused by a sub¬ 
set of carbon atoms in the buffer that locally bond to the 
SiC surface. 
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